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A general review is presented on the problem of non perturbative computation of the K 
amplitude. 
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1. INTRODUCTION 

Kaon Physics is a very complicated blend of 
Ultraviolet and Infrared effects. 

The physical problem is the large enhancement 
shown in if-decays: 
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Due to the difficulty of putting the Standard 
Model on the lattice for technical {Mi\r,z » ^) 
and theoretical reasons (the presence of the Wil- 
son term breaks the chiral gauge invariance) we 
must integrate analytically the heavy degrees of 
freedom {W s, Z s). 

Using second order Weak Interaction pertur- 
bation theory, we get for the effective low energy 
hamiltonian density of non-leptonic decays: 



H 



AS=1 



giy I d^xD{x;Mw)T[J^{x)J+^m + 



where D{x\ Mw) denotes the free W propagator: 
dt^p exp(ipx) 



D{x] Mw) 
and: 



(27r)4(p2 + M^) 



(3) 



Om = ("Is + md)sd + {rns - md)s-f^d + h.c. (4) 

Cm is fixed, in Perturbation Theory, so that quark 
masses are not modified by Weak Interactions. 



The term proportional to Om does not con- 
tribute to the physical transition amplitude (a 
matrix element with zero four momentum trans- 
fer, — 0), being a total four-divergence. 

Since Mw >> ^qcd we can use the Operator 
Product Expansion, to get: 

HlLi = [c+(a.)OW(m) +C-(m)0(-Ha*)] (5) 

where A„ = VudV*^ and: 

= [(57^rf)(«7» ± {S7f:u){u^^d)]-[u ^ c] (6) 

In Eq. Q the subtraction point fi is chosen so that 
/i > rric and we are in the favorable situation of 
propagating charm with the consequent GIM can- 
cellation 0^~^ is an operator contributing to 
pure AI — i transitions, while 0^+^ contributes 
both to AI = i and AI — | transitions. The C±, 
computed in Perturbation Theory, show a slight 
enhancement S : 



C-{p^2 Gev) 



C+(Ai ~ 2 Gev) 



(7) 



The rest of the enhancement (w 10) should be 
provided by the matrix elements of O^^-* and is a 
non perturbative, infrared effect. 

In this talk I will be concerned with the Wilson- 
like fermion formulation. Several interesting pa- 
pers deal with problems related to non-octet K- 
transitions, both within the Wilson and staggered 
fermion approaches Q. 

The difficulty of the problem consists, first of 
all, in giving the correct ultraviolet convergent 
dcffiiition of the operators O^^^ . It is well known 
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that, in order to construct a finite composite op- 
erator of dimension 6, Og{^), we must mix the 
original bare operator, Og(a), with bare operators 
of equal (Og*''(a)) or smaller (03(a)) dimension, 
in general with different naive chiralitiesQ. 
In the present case we have, schematically: 

06(m) = Z{fia) X (8) 
X [06(a) + ^ a{go)oi'\a) + ^^O^ia)] 

i 

Being QCD asymptotically free, one could 
think that the mixing coefficients in Eq.[|| could 
be reliably computed in Perturbation Theory. 

Numerical attempts (and theoretical consider- 
ations) show that this is not necessarily the case. 
We are therefore led to employ general non per- 
turbative techniques. On very general grounds, 
these techniques are based on the systematic ex- 
ploitation of (continuum) symmetries. In fact it 
turns out that the O'^'^^'s have definite transfor- 
mation properties under the SU{'i)®SU{i) chiral 
group (and some discrete symmetries). 

This fact suggests the use of Chiral Ward Iden- 
tities (or equivalent methods) to classify the com- 
posite operators 

The Chiral Ward Identities, however, do not 
fix in a completely unambiguous way the O^^^ s. 
Some ambiguities are left, outside of the chiral 
limit. These ambiguities turn out to be irrele- 
vant in cases in which Current Algebra may be 
applied. They, however, limit the application of 
the method to light meson decays. In order to 
treat other interesting physical processes, involv- 
ing e.g. B-mesons, other methods must be used, 
as discussed later. 

2. ULTRAVIOLET PROBLEMS 

2.1. The Chiral Limit 

As already stated, in order to construct an 
operator with the correct (continuum) chirality, 

we must mix it with operators Oa"^ of equal or 
smaller dimension, with different naive chiralities: 



0„ 



,.n(i 
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A general way to see why Eq.Q is correct, is 
to think of the lattice discretization of QCD as a 



chiral violating order a perturbation of the con- 
tinuum field theory, described by an action: 



Sqcd — S 

cont 

j4 



d yq{y)mq{y)+a / d yW{y) 



(10) 



where Scant describes the continuum theory, in 
the chiral limit, and W{y) is a chiral violating, 
dimension 5 operator whose origin can be traced 
back to the presence, in the Lattice Action, of 
the Wilson term. The existence of W{y) forces 
us, already in the chiral limit, m = 0, to perform 
a complicated subtraction procedure in order to 
define finite ultraviolet composite operators with 
good chiral transformation properties. 
Formally we may write: 

(A|O„(0)|i?)(,,,)«(A| O«(0)|i?)(_,) + 
+ ^ ^ (A| TO„(0)( / dyWiy)r \B) (11) 

where Oq.(0) is, by definition, a finite operator 
which satisfies the chiral continuum (integrated) 
Ward Identities: 



hm / dz(A|r(^„p[g-(z)mg(z)]O„(0))|i?)(_,, 



{A\Sopdc.{0)\B) 







(12) 



where Sop denotes the operator chiral variation. 
The Ward Identity Eq.||l^ is written in a very 
general form and can be applied both to the 
case in which chiral symmetry is restored in the 
massless limit, and to the case, realized in QCD, 
where chiral symmetry is spontaneously broken 
as m — > 0. In this case the mass insertion term in 
Eq. ||l2| does not vanish in the chiral limit because 
of the presence of the pion pole. It is to be noted 
that no prescription is needed for the (generally 
divergent) insertion of Sop[q{z)mq{z)] and Oq,(0), 
since we are considering the limit m ^ and the 
divergent part does not contain the pion pole. 

The a dependence cannot be directly read from 
Eq. ||ll| because of the presence of multiple in- 
sertions of composite operators Oq(0) and W{y) 
which are, in general, divergent. 

For example for n = 1 we have: 

(^|O.(0)|i?)(,,,)«(A| O„(0)|i?)(_,) + 
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+a j d2/(A|T(O„(0)M^(y))|i?)(_,) (13) 

The (continuum) Operator Product Expansion 
implies (up to log corrections): 

1 



T{d^{0)W{y)) « V 



(5+do-<i (i)) 



y 



O«(0) (14) 



Because of the multiplying a factor, in Eq.[ll3| 
only the short distance behavior contributes at 
leading order in a: 



(^1 0.(0) IB) 
«(A|O„(0)|B)(_,) 



(15) 



1 



{A\O^H0)\B) 



[cont] 



The oi*^s in Eqs.||lJ],|l5j are operators which 
generally belong to chiral representations differ- 
ent from that of Oa (apart from trivial multiplica- 
tive renormalisation). In the Chiral Limit the 
Ward Identities completely determine the struc- 
ture of the composite operator Oa- 

In fact Eq.p^ implies the mixing structure 
which has the form indicated in Eq.|^: 

(^1 0.(0) « {A\ O„(0) ^16) 

"E;k;^(^|o«(o) IB) + 

Oa satisfies the Ward Identities up to 0(a). 

We also see that the leading divergent part of 
the mixing, -(35-ra— is determined by pertur- 
bation theory (short distance). 

2.2. Softly Broken Chiral Symmetry 

When the explicit chiral symmetry breaking is 
switched on (to 7^ 0), other divergences (and am- 
biguities) may arise. 

This happens also in the continuum. 

In order to discuss the case to ^ it is more 
suitable to think of a generalized mass term of 
the form qmMijqLj + h.c. (in the end we shall, of 
course, choose M as an hcrmitian diagonal ma- 
trix) . 

If we transform: 



QLi LiklLk 



and, simultaneously: 
M RML+ 

then the mass insertion remains invariant. 
Up to first order in M, we have: 



(18) 



{A\Oa{Q)\B)^^,)^ {A\Oa{Q)\B) 

(M=0) 

dy {A\ Tdam{y)Mq{y) |i?)(M=o) (19) 



As usual the double insertion of composite op- 
erators gives rise to a divergent contribution de- 
termined by the operator product expansion: 



T{Oa{0) q{y)Mq{y)) 
1 



E 



y(?, + dfy-dA-) 



y«0 

(MA)W(O) 



(20) 



(17) 



where the operators ^^(0) have chiralities differ- 
ent from that of 0.(0). {MA)^a denotes a suit- 
able combination of the mass parameters M and 
the operators ^^(0). 

Under a simultaneous chiral rotation, Eq. JTzt, 
on both sides of Eqj2G| and the corresponding 
transformation Eq. |18[, we see that the operator 
combination {MA) a transforms according the 
same chiral representation as 0.(0), since qMq 
is invariant 1^: 

5{MAt^ = {5M A)« + (M 5„pA)« « 5opda{2\) 
We can now define the finite operator: 

ba{Q) ^ 0.(0) - ^ -^^^-^--^(MA)«(0)(22) 

i 

and study its Ward Identities. 

We start with the (divergent, but regularised) 
Ward Identity: 

dy {A\ T{M,,6opmy)qj{y)]dam \B) + 
{A\SopdaiO)\B) = 

dy (A\ Tiq{y)iSM)q{y)dam \B) + 

+ {A\6opda{0)\B) =0 (23) 

Divergent contributions are present in Eq. psf . 
They arise from: 
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1) the double operator insertion in the T- 
product 

2) since 5op denotes the operator chiral varia- 
tion, Oq(0) transforms as a superposition of the 
representations of Oa (0) and of A\ (0) : 

5„pda(0) = ,5opO„(0)+ (24) 

i 

SO that the single insertion of 5opOa{0) is also 
divergent. 

In order to make everything finite we may rede- 
fine a T*-product according to the prescription: 

T*((5„(0) q{v)Mq{y)) = (25) 
= T{b^{0)q{y)Mq{y)) + 

i 

In this way we rewrite the Ward Identity 
Eq.H as: 

- j dy {A\T*q{y){5M)q{y)~d^{Q)\B) + 

+ {A\5opb^{Q)\B)+ (26) 

i 

From: 

- {SM A)« = (M ^„pA)« - 5(MA)« (27) 
we then get: 

- j dy {A\T*q{y){5M)q{y)b^{Q)\B) + 

+ {A\5opbo.{0)\B)+ (28) 
+ E ,fe-l,-i) {A\{M5.pAi^mB) + 

i 
i 

Eq.[^ is easily seen to be equivalent to: 

J dy {A\ T*{M,,Sop{q^{y)qj{y)) O^O)) \B) + 



+ {A\Sd^iO)\B) = 

^-Jdy {A\ T*{q{y){SM)q{y) (5„(0)) \B) + 

+ {A\6 dc,{0)\B) =0 (29) 

So that, in the end, with a joint redefinition of the 
operator and the T-product containing the mass 
insertion, we managed to have a renormalised 
Ward Identity of the canonical form. It is now 
clear why the Ward Identity does not uniquely 

define the finite operator Oa{0), outside the chi- 
ral limit. In fact we may change the definition of 
the T*-product and accordingly change the def- 
inition of the operator Oq(0) by a finite mixing 
with (MA)a'' , leaving the form of the Ward Iden- 
tity unaltered. 

This ambiguity is harmless if we are in a posi- 
tion to apply the low-energy theorems of Current 
Algebra (small quark masses) or when the phys- 
ical matrix elements of {MA)a'' are zero (e.g. if 
they are proportional to four divergences of cur- 
rent), but it can represent a serious obstruction 
when heavy hadrons are present. 

We conclude this section by remarking that the 
situation presented here is a bit simplified. In 
fact we should also consider multiple insertions of 
q{x)Mq{x) and W{y), which, in general, generate 
further distortions. Such terms will be included 
in the following considerations. 

2.3. Non Perturbative Renormalisation 

A different, but equivalent construction of 
composite operators (in agreement with the 
Ward Identities for physical applications) can be 
achieved through the so called Non Perturbative 
Renormalisation 01 . 

This method proceeds in complete parallelism 
with the corresponding procedure used in pertur- 
bation theory. 

In fact non perturbative renormalisation re- 
quires a non-perturbative gauge fixing procedure, 
usually in the form of a suitable discretization of 
the Landau gauge d^Gp = [|. 

^ Gribov copies should not present a serious problem in the 
present case, because gauge fixing is only used in the com- 
putation of short distance quantities which are probably 
immune from their presence. 
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The basic idea is to consider the insertion of 
the renormaUsed operator Oq — Oq + X^^iOq^ 

i 

in Green's functions containing elementary quark 
(and gluon) fields. 

Symbolically we write: 



(30) 



G.(p)=F.T. (O^i0)qix)qiy) 



where F.T. denotes the appropriate Fourier 
Transform of the external legs. 

We expect that, for large (continuum) virtual- 
ities ^ » p ~ fi » Aqcd, the chiral violating 
form factors of (p) (which come from physical 
soft mass breaking and/or spontaneous symmetry 
breaking) should be suppressed by inverse powers 
of 11. 

As an example we could consider the correlator 
F.T. {s{x) SRdL{0) d{y)y In this case the condi- 
tions are: 



F.T. {sl{x) SRdLiO) dLMv)) 
F.T. (slM^) SRdLiO) dniy)) 





(31) 



This behavior is a direct consequence of the (in- 
tegrated) Chiral Ward Identity. I will discuss, for 
simplicity, the case in which mixing with opera- 
tors of lower dimension do not occur. This case is 
directly relevant for the study of the transitions 
K+ TT+TfO and K° ^ K° . 

The integrated Ward Identity has the form: 



F.T. I dz(5oM^Wq{z)]d^{Q)q{x)q{y) 



+F.T. {^S Oo.{0)qix)qiy)) + 
+F.T. (^da^iO) Sqix)qiy)) + 
+F.T. (dc.{0)q{x) Sq{y))^0 



(32) 



Because of the explicit M factor, the l.h.s. 
of Eq. [^2| has one operator dimension less than 
the individual terms appearing on the r.h.s. so 
that, at large virtualities, it vanishes one power 
faster. For asymptotically large n we then get 
from Eq.il: 



F.T. (6 o^Mqi^My) 



+F.T. {O^{0)Sq(x)q(y)j^ + 



-F.T. 
E F.T. 



Oa{0)qix)Sq{y)) = 



6 [da.iO)qix)qiy)] 







(33) 



which is equivalent to the Wigner-Eckart Theo- 
rem (at large virtualities): 



F.T. 



r,d^{0)q{x)q{y) |0) = 



(34) 



Thus, adjusting the mixing coefhcients c; so 
that the chiral violating form factors vanish 
asymptotically, we get an answer equivalent to 
the application of the (softly broken) chiral Ward 
Identities. 

Apart from its simplicity, a further advantage 
of the non perturbative renormalisation approach 
is that it can also determine the absolute nor- 
malization of Oq(0), corresponding to the one 
adopted in perturbation theory, by imposing the 
same kind of conditions as, for example: 



1 



(35) 



F.T. l^qqqqOa 

3. INFRARED PROBLEMS 

In order to compute the K — > tttt width 
we have to evaluate the matrix element 
(out) {''^{p)''^{~P)\ Hw \K) with two interacting 
hadrons in the final state. 

This is not easy to do in the euclidean region 
In fact it can be shownQ that the usual 
strategy of taking large time limits of euclidean 
correlators does not give direct information on 
(out){''^{p)''^{~p)\Hw\K)^ but is contaminated 
by final state interaction effects: 



TTp{h)n^p{t2)Hw{0)K{tK) 



(36) 



tK — > — 00 Z 
tl>>t2»0 

{(out) {7T{p)n{-p) \ Hw \K) + 
+(^n){^T{p)^T{-p)\Hw\K) + 

4_p'^{E^-m„)t2 nj(o.s.) \ 

where m'^^t^'I, tttt is the (tttt — > tttt ) off-shell 

scattering amplitude representing the final state 
interaction. 
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If j3 = we have: 



(37) 



x(l 



tl>>t2>>0 

mKtK-m^t^-m^t2 (7r(0)7r(0) | if^/ \K) X 

c 



^ TTTT — > TTTT ) 



Eq. shows that only for tt s at rest (and 
weakly interacting) it is possible to extract a 
meaningful matrix element. 

For if —> TTTT this could work because chiral 
symmetry implies a small final state interaction 
(Adler zeros). It is however a subject which cer- 
tainly deserves further study. 

We conclude remarking that the euclidean 
Green's function for the A/ = i transition, de- 
fined in Eq. contains, in general, a discon- 
nected contribution 

(^o(ii)^o(t2)) {HwiO)K{tK)) (38) 

which must, of course, be subtracted in order to 
get the physical result. 

4. HOW TO COMPUTE K ^ mr 

This section is based on the considerations ex- 
posed in Ref.|^. 

We will discuss several possible ways to com- 
pute K tttt: 

• The Direct Method (and variants thereof) 
(Parity Violating) 

• The Parity Conserving Method 

• First Principles 

4.1. Direct Methods 

These methods use the strategy first envisaged 
in Ref. [0. The relevant matrix element is eval- 
uated for TT s at rest, which minimizes the final 
state interactions. 

Flavor and CPS-symmetry (CP (g) (s ^ d)) 
severely constrain the form of the subtractions 
needed to define the Parity Violating renor- 
malised Weak Hamiltonian density: 



O 



[P.V.) 



^(P.V.)(Ma) X 



(39) 



C 

^('-'(p.V.)(«) + ("^c - m„)^— Op(a)) 



where: 

Op = {nis - nidjs-i^d (40) 

If we choose a world in which = rrid we 
see from Eq. that the power divergent mixing 
with Op is eliminated. 

In this situation no subtractions are needed. 
The absolute normalization can be found both 
perturbatively or non perturbatively, e.g. using 
external quark states, as in Eq.[^p|. 

We can therefore compute the zero three- 
momentum transfer matrix element: 

A^^lm. = (7r(0)7r(0)| 0^^\^,) |if(0)) (41) 

ms—nid 

and reach the real world through chiral perturba- 
tion theory: 



Ai±) 
phys 



m 



K. 



TT. phys 



(42) 



Another possible way to eliminate the Op sub- 
traction in Eq. |]39| is to use a zero four- momentum 
transfer matrix element. The basic idea is to work 
with non perturbatively 0{a) improved fermions 
(i.e. an S-W improved action with a coefficient 
csw determined as in Ref . [Tl|] ) , and choose the 
quark masses so that rriK — 2m.^. 

In the continuum, since we are working at zero 
four-momentum transfer, we would have: 



{h\ Op \h') cx (tos -I- md) {h\ s^^d \h') 



(43) 



The computation performed on the non- 
improved lattice would give: 



{h\Op\h') 



cx (tos -I- rud) 



{h\sj5d\h') 



{h\d^{sjp^5d)\h') {h\XA\h') 



= (44) 
0(1) 

a a 
and therefore the unsubtracted result would be 
affected by an error 0(1). 

Working with improved axial current and pseu- 
doscalar density, we would have: 



{h\ Op \h'[ 



cx (TOs 



md) 



{h\s^,d\h'\ 



{h\d^isjpj5d)\h') , 0(a2 



a 



^The disconnected contribution, Eq. | 
symmetry. 



0(a) (45) 



j8 , vanishes by CPS 
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In this case we may then compute directly: 



'^(pc)(a^) ^ rather complicated ultraviolet 



(46) 



without worrying about the subtraction. 

The connection with the real world is again ob- 
tained through chiral perturbation theory: 



A, 



(±) 

phys 



m 



phys 



mi 



phys ^(±) 
2 ^rriK^ 



(47) 

'K, lat "^TT, lat 

In order to eliminate the disconnected contri- 
bution, Eq.[Q, we may work with finite, inter- 
mediate renormalised operators obtained by ad- 
justing Cp^^ so that: 



{o\^[p,y,^{^i)\K) = o 



(48) 



This subtraction does not change the decay rate 
because it amounts to change the operator by a 
four-divergence and reduces numerical errors. 

4.2. Real World 

Still another possibility, based on the proposal 
of Rcf. to deal with final state interactions, is: 

1) parametrize the real world K — > tttt ampli- 
tude in the Minkowski region, in a way apt to 
analytic euclidean continuation; 

2) try to fix the parameters by a numerical 
study of the lattice euclidean correlation function. 

This method would require to deal with finite 
operator insertions (also for Ap ^ 0) to get a 
smoother extrapolation. This is accomplished 
through the same pre-subtraction described in 
Eq. ji^ needed also to eliminate disconnected con- 
tribution, Eq.|Q. Since no chiral extrapolation 
is needed, this method could be suitable to the 
study of B-decays. 

4.3. Parity Conserving Method 

This was the method originally proposed in 
Ref. 1^. It is based on the fact that chiral per- 
turbation theory relates the physical matrix ele- 
ment: 



(7r+^-|0gV.)(M)|^")=*7 



ml^ - mz 



(49) 



to the matrix element of the parity conserving 
operator: 



structure: 

x(±) 



O 



(P.C.) 



(51) 



i=l 



(±) 



+ {mc - TO„)C^'''^Oo-(a) + {rric - m„)^^Os(a) 



where: 

Off = goims - rnd)s(7^yG^i,d 

and: 

Os ^ sd 



(52) 



(53) 



however, in this case we only have to compute two 
particle matrix elements (which are easier) and, 
therefore, we do not have to worry about final 
state interaction effects. 

4.4. First Principles 

The distortion of the composite operators in 
the effective Weak Hamiltonian density is due to 
the fact that, in its computation, we insist to in- 
tegrate up to zero (lattice) distance. 

If we could integrate up to a finite, physical, 
distance and perform analytically the final contin- 
uum integration the problem could be overcome. 

In order to accomplish this, we recall that 
the Weak Hamiltonian density Eq. Q is defined 
through a short distance expansion. 

We can therefore directly get the continuum 
matrix elements {h\ 0^^\^ji) \h') by studying, nu- 
merically on the lattice, in the region a « 
|x| << Aqp^, the quantity: 



\[{h\T[J'i^{x)Jt^m\h')+^ 
= 5]c.(x;a.) {h\O^Ht^)\h') 



(54) 



where the coefficients Ci{x] fi) are related to 
the Ci{j^) appearing in the Weak Interaction 
hamiltonian density by: 



{n+ip)\0\p^^{p) \K^iq)) = -5(±)!^+^(±Vg(50) a{^)M^''^ = / dxD{x;Mw)c.{x;pi) (55) 
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The Ci{-^) s, and therefore the Ci{x] n) s, are 
rehably computable within perturbation theory. 
If we use in Eq.[Q the 0{a) improved weak cur- 
rents, then the whole computation will be auto- 
matically 0{a) improved. 

In the case oi K tttt, the operators contribut- 
ing in Eq.|^^ are O^^'' and (m^ — m^)Op, in the 
parity violating case, or (m^ — m^)(ms -I- md)Os, 
in the parity conserving case. 

The perturbative expression of the Ci{x; fi) is: 

Ci{x;fi)cc (56) 



where the anomalous dimensions 7q*'' are given 
by: 



7 



(+) 



4 7 



(-) 



j(P'') = 16 (57) 



As with the previous method, no chiral expan- 
sion is necessary and therefore it could be applied 
also to B-decays. 

4.5. Top Quark and CP Violation 

The difficulty here is that the top quark is es- 
sential and the GIM mechanism is not operative. 
Even in this case, in analogy with the treatment 
described in the previous section, it is possible to 
work in a world with a light fake top quark in 
order to fix the matrix elements of the relevant 
operators. 

In fact we can consider a world in which the top 
quark has a fictitious mass rht such that ^/^ >> 
ifit >> rric » ^QCD so that we can choose the 
subtraction point both in the region ^/^» ^ » 
rht » ™c >> ^QCD where GIM is operative, 
and >> rht >> M >> ''^c >> ^qcd, which 
mimics the real world. 

Identifying the two matrix elements (appropri- 
ately evoluted in /x and studied as functions of 
rht) we can then extrapolate the transition ma- 
trix elements to the real world 
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